Bose-Einstein Condensation in solid ^He 
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We present neutron scattering measurements of the atomic momentum distribution, n(k), in solid 
helium under a pressure p = 41 bars and at temperatures between 80 mK and 500 mK. The aim is 
to determine whether there is Bose-Einstein condensation (BEC) below the critical temperature, Tc 
= 200 mK where a superfluid density has been observed. Assuming BEC appears as a macroscopic 
occupation of the fc = state below Tc, we find a condensate fraction of no = (-0.10 ±1.20)% at 
T = 80 mK and no = (0.08±0.78)% at T =120 mK, consistent with zero. The shape of n(k) also 
does not change on crossing Tc within measurement precision. 



In 2004, Kim and Chan reported the spectacular 
observation of a superfluid density in solid helium below 
a critical temperature, Tc . The superfluid density, ps , is 
observed as a non-classical rotational inertia (NCRI) in 
a torsional oscillator (TO) containing solid helium. The 
percent of solid that has a NCRI and is decoupled from 
the rest of the solid was ps{T) ~ 1.5 % at temperature 
T = 50 mK in commercial grade purity ^He which con- 
tains typically 0.3 ppm of '^He, where Tc — 200 mK. All 
other impurities are frozen out. A ps was observed in 
both bulk solid helium 2] and in solid confined in porous 
media (Vycor) [l]. The magnitude of ps varies somewhat 
from solid sample to solid sample A superfluid den- 
sity was observed in solids at pressures between p 25 
bars near the melting line and 135 bars with ps taking 
its maximum value of 1.5 % at p 50 bars. 

This remarkable result has been confirmed in indepen- 
dent TO measurements 0, S @| ■ Rittner and Reppy Q 
find that ps can be significantly reduced by annealing the 
solids near their melting temperatures with ps reduced to 
zero in some cases. Similarly, Shirahama and coworkers 
[3| report a reduction in ps of up to 50 % by anneal- 
ing. Macroscopic superflow was not observed in helium 
in Vycor Q and bulk helium However, Sasaki et al. 
have observed macroscopic superflow in those solids 
which contain grain boundaries that extend across the 
solid. This unexpected result suggests that there is in- 
deed superflow and that it is along or associated with 
rain boundaries. Superflow related to grain boundaries 
the variation of ps from sample to sample Q and the 
reduction of ps by annealing [1, IH] suggest that a super- 
fluid density may be associated with macroscopic defects 
that extend across or whose impact extends across the 
whole solid. 

In liquid helium, Bose-Einstein condensation (BEC) 
and superfluidity are observed together. Indeed super- 
flow can be shown to follow from BEC 0, [l]| . It can 
also be shown to arise from long range atomic exchanges 



[l^ . In the liquid, both in bulk [l3| and in Vycor [l^, 
BEC is observed as a macroscopic occupation of the k = 
state in 7i(k) , as expected for a translationally invariant 
system. Solid helium is also translationally invariant, as 
is a gas of vacancies in the solid. Thus, for certain mod- 
els of superflow we anticipate that BEC will appear as a 
macroscopic occupation of fc = in n(k) . In this con- 
text we look for an enhancement of n(k) at ~ below 
Tc ■ We also look for a change in shape of n(k) below 
Tc ■ The central result is that we observe no increase 
in n(k) at fc ~ nor any change in shape of n(k) as the 
temperature is lowered below Tc . 

In 1969, Andreev and Lifshitz [16] proposed that he- 
lium could be a supersolid if the solid contained vacant 
sites in the ground state. Essentially, the ground state, 
zero point vacancies could form a Bose gas at T ~ K 
in which the Bose-Einstein condensate fraction, no , and 
ps are both approximately 100 %. While thermally ac- 
tivated vacancies have been observed [l3|, ground state 
vacancies have not. Chester [l^ proposed that super- 
flow in solid helium may be possible because the solid 
is well described by fluid like wave functions that sup- 
port superflow. Leggett (19| examined superflow via long 
range exchanges of atoms within a perfect helium solid 
and found superflow possible but that ps would be small, 
PS - 0.01 %. 

Recent accurate path integral Monte Carlo calculations 
^2^ find that ps arising from long range atomic exchanges 
in a perfect crystal will be unobservably small. Similarly, 
BEC in bulk perfect crystals is predicted to be very small 
22l | . Zero point vacancies in the ground state are pre- 
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dieted to be unstable [2l| . Essentially, vacancies migrate 
to a surface or coalese to create a surface and leave the 
crystal. However, if solid helium is held in an amorphous 
rather the equilibrium crystal state, both ps and no take 
significant values [2l| . With these predictions, it is inter- 
esting to search for BEC where ps is observed. 

The solid ^He investigated in this experiment was 
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FIG. 1: Observed J{Q,y) of solid ^He at V = 20.01 cm3/mol 
(open circles) folded with the instrument resolution at momentum 
transfer Q = 26 and temperatures T = 120 (top) and 300 

mK (bottom). The solid lines are fits of the Convolution Approach 
(CA) with Final State (FS) function taken from liquid *He (Ref. 
|l3|). shape of n(k) from Ref. [24| and width 02 parameter as the 
single free fitting parameter. The data peaks below y = because 
of FS eff'ects. The dotted line is the MARI instrument resolution 
function. 



grown using the blocked capillary method. Commercial 
grade purity "^He (0.3 ppm '^He) was introduced into a 
cylindrical Al sample cell of 20 mm diameter and 62 mm 
height at a temperature T ~ 3 K to a pressure of p ~ 70 
bars. At constant p, the temperature was reduced using 
an Oxford Instruments Kelvinox VT dilution refridgera- 
tor until solid formed in the capillary and blocked the cell. 
The block was observed at filling p — (69.8 ±0.2) bars 
and T — (2.79 ±0.02) K, which corresponds to liquid on 
the melting line at a molar volume Vm — (20.01 ±0.02) 
cm'^/mol 23]. The blocked cell was further cooled and 
neutron inelastic scattering measurements at high mo- 
mentum transfer were made in the solid hep phase at 80 
mK, 120 mK, 300 mK and 500 mK on the MARI spec- 
trometer at the ISIS neutron facility. 

Single atom dynamics is observed in the dynamic stru- 
ture factor S{Q,lu) at high momentum transfer, HQ. 
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FIG. 2: The mean square atomic momentum, 02 = {kg) along an 
axis and the atomic kinetic energy K = (3A/2)ci2 where A = 12.12 
K of solid helium at V = 20.01 cm'^/mole versus temperature. 

Specifically, at HQ — > oo, where the Impulse Approxi- 
mation (lA) is valid, the energy transfer huj in S{Q,uj) is 
Doppler broadened by the atomic momentum distribu- 
tion n(k) . In this limit, it is convenient to express lu 
in terms of the '2/ scaling' wave vector variable, y = 
(uj — ujpi)/vR where (jJr = fiQ"^ /^m and vr = hQ/m, 
and to present the neutron inelastic data as J{Q, y) = 
vrS{Q,uj). Including final state (FS) effects which are 
small but not negligible at the Q values investigated here, 

J(Q, y)^ J dy'R{Q, y - y')JiAiy') (1) 

where R{Q, y) is the FS broadening function and 

JiA{y) = j d\^n{\^)5{kQ - y) = nQ{y) (2) 

is the lA to J{Q,y). Specifically, J/^(?/)is n(k) projected 
alon g Q denoted the longitudinal momentum distribu- 
tion [25|. 

Fig. [1] shows the observed J{Q, y) at wavevectors Q — 
26 A-i and temperatures T = 120 mK and T = 300 
mK. The observed J{Q, y) includes the MARI instrument 
resolution function which is shown separately as a dotted 
line in Fig. [TJ The solid line is a fit of a model to the 
data as described below. We were able to determine at 
most two free parameters in model fits to the data. 

To obtain a condensate fraction, we assumed a model 
n(k) of the form, 

n(k) =noJ(k) + (l-no)«*(k), (3) 

where n* (k) is the momentum distribution of the atoms 
above the condensate in the fc > states. To proceed, we 
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FIG. 4: Fits of the Additive Approach (AA ) which includes FS 
effects and a simple Gaussian to the observed J{Q,y) ■ 



FIG. 3: Condensate fractions obtained by fitting the Convolution 
Approach (CA) at T = 80 mK and 120 mK with the width of 
n(k) (02 ) held fixed at its values shown in Fig[2] We find no = 
(-0.01 ±1.20)% at T = 80 mK and no = (0.08±0.78)%. 



assume (1) that the shape of n(k)is the same as observed 
previously in soUd hehum at a somewhat lower pressure 
[ii] and (2) that the FS function R{Q, y) in Eq.(l) is the 
same as observed [l3l| in liquid helium. The free parame- 
ters (two) in the model are then uq and the width, 02, of 
the Gaussian component of n* (k) . A condensate compo- 
nent appears as an additional intensity unbroadcncd by 
n*(k) in JiAiu) at y = 0. 

We first fit the model n(k) to the data at 500, 300, 
120 and 80 mK assuming ng = at all temperatures. 
The no is expected to be zero at 300 and 500 niK. The 
resulting values of a2 and the associated atomic kinetic 
energy {K = (3nV2m)a2) are shown in Fig. [2l The 
resulting 02 and K decrease somewhat with decreasing 
temperature. In a recent measurement, Adams et al. |14l | 
find K independent of temperature within precision. 

If superflow is associated with defects in the solid such 
as vacancies, we anticipate that BEG is similarly asso- 
ciated with these defects, perhaps in a vacancy gas. In 
this event, the atomic kinetic energy, K, of the majority 
of the atoms may be largely unaffected by the BEG in 
the defects. To obtain no at 80 mK and 120 mK within 
this picture we kept a2 fixed at the values obtained above 
and refitted model ^ with no as a free parameter. The 
fitted values of no are shown in Fig. [3l The variation of 
no with Q reflects the statistical error in no . The mean 
values are no (-0.1 ± 1.2) % at 80 mK and no = (0.08 
± 0.78) % at 120 mK. 

If, in contrast, the superflow and BEG lie within the 
bulk of the solid, we anticipate, as in liquid ^He, that 
the atomic kinetic energy, K, will decrease below Tc as a 
result of BEG. The observed decrease in K below Tr has 



been used to estimate no in liquid *He 2^ 27 1. To ad- 
dress this case, we kept a2 (kinetic energy arising from 
n* (k) , constant at the value obtained at 300 mK and re- 
fitted model ([3]) at 80 mK and 120 mK to obtain no . 
The resulting values are no = (0.8 ± 1.2) % at 80 mK 
and no = (0.76 ± 0.77) % at 120 mK. This method as- 
sumes that all the drop in K below Tc = 200 mK arises 
from BEG. To normalize these no values for changes in 
K from other sources, we also determined n-o at 300 mK. 
That is, we kept the a2 fixed at its 500 mK value and 
refitted the model to determine uq at 300 mK giving 
no = (0.63 ± 0.77) %. Since no at 300 niK is zero, 
we expect the no values below Tc to be overestimated 
by approximately 0.6 %. Normalizing for this overesti- 
mate, we arrive at no — (0.2 ± 1.2) % at 80 mK and 
no ~ (0.1 ± 0.8) % at 120 mK. If we try to determine 
both a2 and no simultaneously, we get essentially the 
same values as before, with only larger error bars; e.g. 
no = (0.74 ± 1.01) % at T = 120 mK. Thus aU methods 
give similar values of no . 

To investigate a possible change in shape of J{Q, y) or 
n(k) on crossing Tc = 200 mK, we fit the additive ap- 
proach (AA) [28[ to the data. In this model to lowest 
order, J{Q,y) is [H,!!! 



Jg{x) 



(4) 



where Jr' = 



■■exp{—^) and x = y/a2^^^- The sec- 
ond term in (4)is the leading FS term and the third term 
is the leading correction to a Gaussian n(k) plus a FS 
term. The fitting parameters are ^2, /J3 = 0,3 /XQ, and Ji^ 
= cti + ai/ {XQY . Previously we found 04 = [l3l. Iisl. l2a| . 
The three remaining parameters are 1x2,713, and 54 . The 
kurtosis of n(k) is 5 = =^ . 
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FIG. 5: Fitting parameters h^^i2 = (AQ)^ 02 and h^^is = 
(AQ)2 as in the AA giving 02 = 1.67 A'^ g^^j ^^/^ = g 21 A""'. 



Fig. m compares AA and simple Gaussian fits to data. 
Tlie AA gives a better fit witfi the peak position lying 
at y < in agreement with the data. We were able to 
determine only two parameters, e.g. h^fi2 — {^Q)^a2 
and h^H3 = (A(5)2 03 where A = h'^/m = 1.0443 meV 
= 12.12 K A^. These parameters are plotted in Fig. [5] 
and show the expected Q dependence. We found again a 
smooth change in 712 with temperature as shown in Fig[2] 
and found 03 /A independent of temperature. We tried 
keeping 7x2 fixed and fitting for /I3 and 754 but we were not 
able to determine CI4 . 6 — or 5 = 0.4 fit equally well. 
Thus we find no change in the shape of J{Q, y) or n(k) at 
Tc within precision. 

The condensate fraction, nn, of a perfect crystal is ex- 



pected to be very small [21|, [22], no < 10~^. In a per- 
fect crystal, BEG requires double occupation of a lattice 
site which has very high energy and is therefore very 
improbable. If there are vacant sites, ng in the crystal 
is dramatically increased, to no — 0.23 % for a vacancy 
concentration — 1 % [1^. If the solid is frozen in a 
non-equilibrium amorphous state, uq ~ 0.5 % [lO]- The 
no in the amorphous state is relatively insensitive to den- 
sity although ps decreases significantly with increasing 
density 30] . 

If the vacancies are treated as an ideal Bose gas in 
which no~ 100% and ps— 100%, a vacancy concentration 
gas cv ~ 1-5 % is needed to reach ps = 1.5 % in the solid. 



This also predicts no ~ 1.5 %. This large value of no and 
this simple model appears to be excluded by our observed 
values, of no= (-0.10 ±1.20)%. In contrast, the nn values 
within bulk solid helium including vacancies [29| or in 
extended amorphous regions [1^ , whether in equilibrium 
or not, are consistent with our observed value. 

In summary, we have determined the BEG condensate 
fraction in commercial grade purity solid helium at pres- 
sure 41 bars and molar volume 20.01 cm^/mole using 
inelastic neutron scattering. We find a condensate frac- 
tion no = ( -0.10 ±1.20)% below = 200 mK consistent 
with zero. We also find no change in the shape of the 
atomic momentum distribution on crossing . 
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